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Quantification of bacteriophages by real-time quantitative PCR (qPCR) is an interesting alternative to the traditional 
plaque assay. Importantly, the method should in principle be able to discriminate between closely related phages that 
are indistinguishable by most other means. Here, a method is presented that employs qPCR to discriminate and quantify 
ten closely related lambdoid phages of Escherichia coli str. K-12. It is shown that (1) treatment of samples with DNase 
efficiently removes non-encapsidated DNA, while the titer of plaque forming units is not affected, (2) individual phage 
types can be accurately quantified in mixed lysates, and (3) the detection limit corresponds to that of a plaque assay. The 
method is used to quantify individual phage types that are released from lysogens that carry up to three different 
prophages. 



Introduction 

The infection of a bacterial cell by multiple, different bacterio- 
phages (phages) constitutes an important part of phage ecology 
that influences phage evolution in many ways.''^ To study inter- 
actions among phages, methods are required that can discriminate 
and quantify different phage types in the same sample. The 
double agar overlay assay,'' although considered to be the gold 
standard for phage quantification, falls short of what is required 
here. Its capability to differentiate phages is limited as it depends 
on differences in phenotypic markers such as plaque morphology 
or host range of the phages to be distinguished.^ 

Real-time quantitative PCR (qPCR) is an interesting alternative 
and has been successfully used to enumerate phages as well as 
other microparasites.^ '" Concerning phages, current applications 
of qPCR fall broadly in two groups, neither of which requires a 
highly specific assay. Either the specificity of the method is not an 
issue as samples are expected to contain only one type of phage, or 
the method has been deliberately adapted to be unspecific such 
that a broader group of phages is detected." '^ Here, a method is 
presented where qPCR is used to achieve the opposite, namely the 
discrimination of closely related phages in mixed lysates. 

To develop and validate this method, a set of ten lambdoid 
phages of Escherichia coli str. K-12 is used. These phages have 
been used previously to study the competitive interactions 
that occur when two prophages are induced in the same lysogenic 
cell.^ Pairs of these phages can be distinguished by differential 



plating on bacteria that exclude one phage by its receptor 
specificity or immunity group. However, this approach fails to 
discriminate most combinations of three or more phages. This can 
be overcome by a qPCR approach where individual phages types 
are specifically targeted. 

The assay is based on an earlier protocol,' a modified version of 
which has been used previously to quantify variants of phage X 
that did not produce plaques or whose plaques were difficult to 
count.'* It consists of two steps, an initial treatment with DNase 
to remove phage DNA that is not enclosed in capsids (which 
would cause an overestimation of the phage titer), and the qPCR 
itself. While the method performs very robustly, its extension to 
achieve the discriminatory power to accurately quantify individual 
phages in mixed samples has never been explored and required 
further experiments for validation. 

Together, these experiments show that the assay detects all 
viable phages in a sample and clearly discriminates individual 
phage types. The good correspondence to the plaque assay 
suggests that the method can be used by digital PCR'^ for 
absolute quantification without an external standard and for the 
quantification of phages which are difficult or impossible to 
quantify otherwise (e.g., because they form no or barely visible 
plaques, or no cells suitable for plating are available).''*"' 

The method is used to quantify individual phages in lysates that 
were obtained from the induction of lysogens containing either 
one, two, or three different lambdoid prophages. Specific pair- 
wise combinations of these phages have previously been shown to 
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strongly interfere during prophage induction, leading to substan- 
tial declines in productivity for both phages involved.^ It is 
therefore of interest to investigate whether this mutual inter- 
ference also affects other prophages that partake in the induction, 
and whether they can mediate the outcome. 

The method presented here is hoped to be a useful tool to 
conduct experiments in the area of phage ecology and evolution. 
Data on viral diversity in various environments accumulate at a 
staggering pace,'^'^° and it is important to test whether concepts 
from ecological and evolutionary theory can be applied to 
understand the processes that give rise to this diversity.^''^'' 
Another field where the application of qPCR is promising is 
phage therapy, the treatment of bacterial infections with phages. 
The ability to discriminate different phages will allow better 
understanding of phage population kinetics when phage cocktails 
instead of single phages are used.^^'^'' 

Materials and Methods 

Materials. Bacteria, phages, and reagents. The phages used were all 
lambdoid phages of Escherichia coli str. K-12 and were described 
elsewhere (Table l).''^^-^'' They all belong to separate immunity 
groups, except for phages mEp234, and mEp332, which share the 
immunity group of phage "k. All phages were kept as prophages in 
E. coli str. K-12 substr. MG1655. 

Bacteria were grown in lysogeny broth (LB; 10 gl^' tryptone, 
5 gl"' yeast extract, 10 gl^' NaCl). For plates, 1.5% agar was 
added, for top agar, 0.7%. Lysates were stored in SMG (0.1 M 
NaCl, 8 mM MgS04, 50 mM TRIS-HCl pH 7.5, 0.01% gelatin) 
at 4°C. Prophages were induced with mitomycin C (Sigma, 
M4287) at a concentration of 5 (J,g/ml. 



DNase treatment was performed with RQl RNase-free DNase 
(Promega, M6101), using lOx reaction buffer (Promega) con- 
taining 400 mM TRIS-HCl pH 8.0, 100 mM MgS04 and 
10 mM CaCl2, and stop solution (Promega), containing 20 mM 
EGTA pH 8.0. 

Primer design. Primer3 software as included in Geneious Pro 
v5.5^'* was used together with sequence data of all phages^^'^' 
(Refardt D., unpublished data) to select primer pairs that 
specifically target selected phages (Table 1). Annealing tempera- 
ture was set to 60°C and product size to -100 bp. Primers were 
ordered from Microsynth (Balgach, Switzerland). 

Instrument settings and data analysis. Amplification reactions 
were performed with a 7500 Fast Real-Time PGR System 
(Applied Biosystems) using default settings including a melt 
curve. Cycle threshold (ct) values, the number of PGR cycles 
required to amplify a target to a detectable level, were calculated 
at a normalized fluorescence emission intensity threshold of 
ARn = 0.1, which corresponded to approximately 5% of the 
maximum intensity. Statistical analyses were performed using R 
statistical software version 2.11.0.^^ 

Procedure. Preparation of phage lysates. Inoculate single 
lysogenic colonies in LB and grow under continued shaking 
overnight at 37°C. Dilute cultures 1/500 into fresh LB and bring 
back to exponential growth by incubation for 2 h under con- 
tinued shaking at 37°G. Induce lysogens by adding 5 (J.g ml"' 
mitomycin C and continue incubation for approximately 4 h at 
37°C. Dilute lysates 1/10 in SMG, add a drop of chloroform, and 
centrifuge at 2700 g for 10 min. Store lysates at 4°C. 

Removal of non-encapsidated DNA. Add 2 \A DNase (2 units) 
and 1.3 |il lOx reaction buffer to 10 |.d lysate and incubate for 
60 min at 37°C. Inactivate DNase by adding 1.5 1^1 stop solution 



Table 1. Primers that were validated in this study and the phages against which they are targeted 



Name 


Sequence 


Target phage 


Target sequence 


Product size 


qPCR_HK022 


5'-TGCCATCGCCATCAAAACAGGT-3' 
5'-TCATCACGGTTCGCGGTGACA-3' 


HK022 


cl 


96 bp 


qPCR_mEp043 


5'-ACCCGCACGAACGTTACCCG-3' 
5'-TCCCACGCGGATGGATGGACA-3' 


mEp043 


cl 


86 bp 


qPCR_mEp213 


5'-TCGCCCTGAACCCTGAGCCA-3' 
5'-GGAGTGAGGCCGTAGAGCCGT-3' 


mEp213 


cl 


103 bp 


qPCR_mEp234 


5'-GCAATTAGTTGGTGCATGCGGCG-3' 
5'-ATCCCCTACCTCAGCGCGGG-3' 


mEp234 


int 


96 bp 


qPCR_mEp235 


5'-TCAGGGGAGCACTGCAAAGCC-3' 
5'-GCAGGCTGCCCTTGGCAAGAT-3' 


mEp235 


cl 


126 bp 


qPCR_mEp235b 


5'-GGAACCCAGCGTGCTGTTGC-3' 
5'-CAGGGCGCCAGCTGTAACGA-3' 


mEp235 


cro 


120 bp 


qPCR_mEpX1 


5'-AAATGGCCCGGCCTGCACrG-3' 
5'-GCACTCGGCGTATCTCCGGC-3' 


mEpXI 


cl 


105 bp 


qPCR_mEpX2 


5'-TGCTATGACCAAGCTGCCGTTGA-3' 
5'-TCAATGACGGCACCGGTTGGG-3' 


mEpX2 


cl 


96 bp 


qPCRJambda 


5'-GCGTTACCGTTTCGCGGTGC-3' 
5'-TCGCAGCATTGCCCGTCAGG-3' 


X, mEp332 


or/51 (phage X) 


120 bp 


qPCR_Phi80 


5'-GCACCCCGCTTGAGAAAGCCA-3' 
5'-CTGGAAGGCTGCCACCrCGC-3' 


O80 


cro 


120 bp 
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and incubating for 30 min at 95°C. 
Heating also opens phage capsids and 
makes DNA accessible for amplification. 
Dilute treated lysates 10-fold in H2O and 
store at -20°C. 

Quantification by qPCR. Amplify 
samples in 10-|ll reaction volumes includ- 
ing 300 nM of each primer, ix SYBR 
Green PGR master mix (Applied 
Biosystems, 4309155), and 1 |il template. 

After overnight incubation of cultures, 
the complete quantification process can 
be done within one day. Multiple samples 
can be efficiently processed in parallel by 
growing cultures in 96-deep-well plates 
(add a glass bead to each well to ensure 
proper mixing of cultures and seal plates 
with breathable film). 

Validation of the protocol. Capsid 
stability. An experiment was conducted 
to test whether treatment of lysates with 
DNase causes a change in the titer of 
plaque forming units. From each phage, 
18-24 separate lysates were obtained and 
quantified using a plaque assay. Lysates 
where then treated with DNase and 
quantified again with a plaque assay 
before the final heat-inactivation step. 

Significant (p < 0.05, Wilcoxon signed- 
rank test) changes in phage titer were 
observed in three of ten phages (Fig. 1). Of 
these, the largest change (+25%) occurred 
in phage mEp332. This change is moderate 
when compared with the inter-assay vari- 
ance, which is typically larger.'" It is 
concluded that the observed effect is not 
relevant for quantification. 

Efficiency of the DNase treatment. The 
amount of phage DNA in lysates that is 
not enclosed in capsids was measured by 
amplifying 3-6 lysates per phage with 
qPGR both before and after DNase treat- 
ment. Gt values were consistently higher in 
treated than in untreated samples (Fig. 2). 
Samples of phage mEp332 showed the 
largest difference (Act = 7.1), which corresponds to a 100-fold 
excess of non-encapsidated phage DNA relative to DNA in 
virions (assuming an amplification efficiency of 95%). 

To estimate DNase efficiency, lysates were kept at 95°G for 
30 min prior to DNase treatment. When these samples were 
quantified with qPGR, their ct values were found to be at 
least 8.4 cycles above those of untreated samples (Fig. 2). 
Thus, DNase removes at least 99.6% of all accessible DNA 
(assuming an amplification efficiency of 95%). Together, the 
results confirm that treatment with DNase efficiently removes 
non-encapsidated DNA from lysates. Using the numbers 
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Figure 1. Ratio of plaque forming units in lysates of different phages after and before treatment 
with DNase. Eighteen to 24 independently obtained lysates were quantified for every phage. 
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Figure 2. Amount of DNA as measured by qPCR in lysates that were either not treated with DNase 
(dark gray bars), treated with DNase (gray bars), or heated to 95°C and then treated with DNase 
(light gray bars). Error bars are standard errors of 3-6 replicates. 
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obtained, the worst possible case is an 
the titer of phage mEp332 by 3% due 
encapsidated DNA. 

Primer specificity. All primers were tested in replicate 
amplification reactions against lysates of all phages and water 
controls. Generally, primers were found to only amplify the 
respective target sequence. After controlling for differences in 
melting temperature (using a difference of > 1°G from that of the 
target sequence as criterion for exclusion), 18 of 196 negative 
controls yielded detectable amplifications that would have been 
accepted as positive results had they occurred in samples with 
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unknown content. Typically, these arose very late (ct > 35, which 
corresponds to less than two copies of the target sequence) and 
inconsistently among replicates. The exception were primers 
qPCR_mEp235 and qPCR_mEp235b, which very weakly cross- 
amplified samples of phages X, mEp332, and mEpXl, and X, 
mEp234, and mEp332, respectively. Ct values of these cross- 
amplifications were all above 30, which corresponds to less than 
10 phages. Therefore, these cross-amplifications cause a negligible 
bias in the quantification of phages. However, detections of 
very low numbers of phages (in particular phage mEp235) in 
samples of unknown content must be interpreted with the 
appropriate caution. 

Amplification efficiency. Lysates from all phages were quantified 
with plaque assays, treated with DNase and then serially diluted 
down to an expected number of 0.01 phages in the amplification 
reaction. Standard curves were constructed by amplifying at 
least four replicates per dilution step (Fig. 3). Amplification 
efficiencies, estimated by linear regression, were found to be 
between 87 and 102% and were estimated using the data that are 
connected by a solid line in Figure 3. This was done to avoid 
overestimating amplification efficiencies, which occurs when 
dilution steps are included that contain replicates where no 
phages were detected. 

Detection limit. If qPCR is congruent with a plaque assay, a 
single phage genome per reaction defines the largest obtainable ct 
value (horizontal dashed lines in Figure 3). When this limit is 
approached, variation among replicates increases as the actual 
number of phage genomes follows a Poisson-distribution. Five 
percent of all samples with an expected number of three phage 
genomes will be empty. Ten percent of all samples with an 
expected number of 0.1 phage genomes will contain a single 
phage genome, the others none. 

Analysis of the amplified dilution series revealed that the 
detection limit was reached when samples contained an expected 
number between 0.1 and 2.8 phages, depending on phages and 
primers used. This suggests that results obtained by qPCR are 
congruent with those of a plaque assay. 

Example: Prophage induction in single, double, and triple 
iysogens. Four triple lysogens (}t/HK022/mEpX2, }^/mEp213/ 
mEp235, ?t/mEpXl/mEp235, }t/mEpXl/mEpX2), as well as 
all possible single and double lysogens of these phages were 
constructed and validated as described previously.'' Of each 
lysogen, seven independent lysates were obtained and treated, and 
the amount of individual phages was quantified by qPCR. To 
obtain standards against which phages can be quantified, one 
lysate from each single lysogen was quantified with a plaque 



assay and used in a dilution series with five steps and two 
replicates each. 

Productivity as measured in single lysogens is interpreted as the 
baseline performance of a phage. ^ All measurements from double 
and triple lysogens were normalized by expressing them relative to 
the baseline performance of the respective phage. Results from 
single and double lysogens were in accordance with an earlier 
study, where measurements were obtained with a plaque assay. ^ 
Together with the small variation among replicates, this indicates 
that quantification with qPCR provides reliable results (Fig. 4). 

Results were analyzed by first grouping them according to each 
of the four triple lysogens and the corresponding double lysogens 
(which corresponds to the four panels in Fig. 4). Data were then 
log-transformed and analyzed with analysis of variance including a 
single factor that separated all nine titer measurements in the 
different lysogens (which corresponds to the nine box-plots in 
each panel of Fig. 4). Titers of the same phage in lysates from 
different lysogens were compared with Tukey's HSD test at a 
95% confidence level. 

The results show that induction of double lysogens that contain 
phage X together with either phage HK022, mEp213, or mEpXl, 
results in lysates where the titer of both phages is strongly reduced. 
Other double lysogens did not exhibit such a strong interference, 
as always one phage kept its titer close to its baseline performance. 

Interestingly, the negative effect of phage X on phages HK022, 
mEp213, and mEpXl, was removed by the addition of either 
phage mEp235 or mEpX2 (Fig. 4). In these triple lysogens, the 
titer of phage X was still strongly reduced, yet the titer of phages 
HK022, mEp213, or mEpXl increased about a hundred-fold to a 
level similar to that of the corresponding double lysogen that does 
not contain phage X. The results show that the interaction of 
multiple prophages can have unexpected outcomes that are not 
predictable based on observations of individual phages. The 
mechanism behind these interactions remains to be explained. 
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Figures (See opposite page). Standard curves for all primers and their respective target phages. The title of every plot gives the names of the primers 
and phage that have been used. The solid line connects data points to which a linear regression model has been fitted to calculate amplification 
efficiency. The regression is then extended by the dotted line down to an expected number of one phage, after which it is continued horizontally to the 
most diluted sample that has been used. Ideally, all data points should follow this line. 
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Figure 4. Quantification of phages that are released upon induction of lysogens that carry multiple prophages. Reproductive success of every phage is 
given relative to that achieved in a single lysogen. Box-plots summarize the results of seven replicates. Colors indicate data from phage X (yellow), 
phages HK022, mEp213, or mEpXI (blue) and phages mEp235 or mEpX2 (brown). 
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